plotG(vx,vy,char,size,color)

9ce0314e-d16f-4b32-ad49-c3d32b86bb1d 10 6 9 decimal

&[DATE] &[TIME] - &[FILENAME]

vx vy char size color plotG n vx length : plot vx 1 el vy 1 el char size color augment : i 2 n range plot plot vx i el vy i el char size color augment stack : for plot 41 line :

Thiele Table

1. Thiele Table ... Thiele(X,Y,Δ) => 'K'

X Y Δ Thiele ε 1018 -^: n X rows : i 1 n range M i 1 el Y i el : for j 2 n range i j n range dM M i j 1 - el M j 1 - j 1 - el - : dM dM abs ε < ε M i j el X i el X j 1 - el - dM / : if : 21 line for for i 1 n range c i 1 el X i el : c i 2 el M i i el : 21 line for Δ 0 ≡ c M if 61 line :

data set XY

data 0.5 - 107.477 0 101.325 0.5 95.461 1 89.876 1.5 84.559 2 79.501 2.5 74.691 3 70.121 3.5 65.78 4 61.66 4.5 57.752 5 54.048 5.5 50.539 6 47.217 6.5 44.075 7 41.105 7.5 38.299 8 35.651 8.5 33.154 9 30.8 9.5 28.584 10 26.499 10.5 24.54 11 22.699 11.5 20.984 12 19.399 12.5 17.933 13 16.579 13.5 15.327 14 14.17 14.5 13.1 15 12.111 15.5 11.197 16 10.352 16.5 9.571 17 8.849 17.5 8.182 18 7.565 18.5 6.994 19 6.467 19.5 5.979 20 5.529 422 mat :

ListCollect(data,select)

data select ListCollect collects index list ct 1 : i select a ct el data i 1 el : b ct el data i 2 el : ct ct 1 + : 31 line for a b augment 41 line :

Cfr(j,K,X,x)

j K X x Cfr Optimiz Symbolic j rows(K)-2 N j 1 + el K j 1 + el x X j 1 + el - K j 2 + el / + : c j el N j 1 + el : a j el K j el x X j el - c j el / + : jj j 1 - (1 range c jj el a jj 1 + el : a jj el K jj el x X jj el - c jj el / + : 21 line for a 1 el 71 line :

Error Stem Plot

Error stem plot for QuickPlot unicolor data Error stemX 1 data rows range : i 1 stemX rows range stem i el data i 1 el data i 2 el data i 1 el 022 mat : for i 1 stem rows 1 - range stem i el stem i el stem i 1 + (1 el 1 el 012 mat stack : for V stem 1 el : j 2 stem rows 1 - range V V stem j el stack : for V V data data rows row stack : 71 line :

Select the support points select 26101418222417 mat :

XY data select ListCollect :

Plot data & support points

X XY 1 col : Y XY 2 col : Kvec X Y 0 Thiele : 31 line

Pts i 1 data rows range p i el . : for data 1 col data 2 col p augment 21 line :

X Kvec 1 col : K Kvec 2 col : 21 line

support X Y + 15 blue plotG :

Support points XY 0 101.325 2 79.501 4 61.66 6 47.217 8 35.651 10 26.499 11 22.699 72 mat

Thiele 'K' convergents Kvec 0 101.325 2 0.091642 - 4 217.335918 - 6 0.210251 8 318.568554 10 0.154456 - 11 136.240642 - 72 mat

Activate Thiele

j Kvec rows 2 - :

j 5

x g j K X x Cfr :

Relative residuals for the "ErrorPlot" Residuals sub data 12412 submatrix : i 1 sub rows range δ i el 1 sub i 2 el sub i 1 el g / - eval : for sub 1 col δ augment 31 line :

Pts x g support 3 1 sys

Relative residuals explorer

Residuals Error Residuals 2 1 sys

Export an easy 'J(x)' fraction ...
evaluate coefficients from symbolic ... 6 decimals adequate. x J 66.316994 3988.071258 x 2.403022 + 1379.737008 x 28.965069 + 539.638243 x 13.217861 + / - / + / - :

Optimiz  ... NONE x g confrac x convert maple 13013516466632178249516061132139237785453196232002317333365704175404275320000000 / 307137309728143861325097941827548993452402936164216938009509465411051459996917701399746693985481551484733870146450997672817982072106872220480000000000 / x 84944583815526841869391889719572753464070919773972898618715991007114136210597927518218738776361387601209690973534907288815991899190147092969855976893470856000308370374608339334526136638013027755523391626664279647464000 / + 279828996928380066397673128112292182096159264977201444161293009928700757832533718034254262449974123233493218084720742295891234371403901775138171777868376202813286407585447740200451522878505658443398089657729075374950471722742425299195485444308160185041366570283421984812716970277681225510341094637080025 / x 17493913050235941722686314608138518624995841499658488142184918324225400912336277573048394754495006085205595578712641393338255512586937421879602488816392807153456832236419098638113092014195262745913741946084447064052603965859535369390037093696084837240877583904627148749650171553874931082715776164435837139367947236822585375287410320996858892589305526339346935826252436659744703214265656533638812023290411986354068910832743247965 / + (9705794607807234932709320075934972990135215167532079781580850606050759002995635316199517741797924455148969259635582436601024445751044494290849836899747331313328326044958532227723205702264344733795380843980355583037850083643799260509662063712346781016578472800000000000000000000001798574274622348199332380121487642613171861394921920173514911924663539175747551789913491984725793834359122264784544568567853808161803546500904379075828671974361955603099561620109262395386961259896716474020668406026383077285036140286918459081276962922348640236162366495054228849 / x 177265966188542847467212659589951892469807868480170359165569183984071922608615686289642834225403598609235538144453425423538449385377264502501341109344767363925973635124639759442948537089827343071386021528799281384878847496968875176289389366512221889189189817969813586907363069543 / + (/ - (/ + (/ -

66.316994 3988.071258 x 2.403022 + 1379.737008 x 28.965069 + (539.638243 x 13.217861 + (/ - (/ + (/ -

χ X Mean : μ Y Mean : n XY rows : 31 line

Pearson’s correlation coefficient corr X i el χ - (Y i el μ - (* (i 1 n sum X i el χ - (2^(i 1 n sum Y i el μ - (2^(i 1 n sum * sqrt / X i el χ - (X i el J eval μ - (* (i 1 n sum X i el χ - (2^(i 1 n sum X i el J eval μ - (2^(i 1 n sum * sqrt / / : 1

X Y 1 Thiele 101.325 000000 79.501 0.091642 - 00000 61.66 0.100845 - 217.335918 - 0000 47.217 0.110889 - 207.823476 - 0.210251 000 35.651 0.121814 - 198.862646 - 0.216529 318.568554 00 26.499 0.133643 - 190.470942 - 0.223339 305.619843 0.154456 - 0 22.699 0.139903 - 186.487519 - 0.226916 300.026627 0.161795 - 136.240642 - 77 mat